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It has been pointed out that the axion electrodynamics exhibits instabihty in the presence of a 
background electric field. We show that the instability leads to a complete screening of an applied 
electric field above a certain critical value and the excess energy is converted into a magnetic field. 
We clarify the physical origin of the screening effect and discuss its possible experimental realization 
in magnetic materials where magnetic fiuctuations play the role of the dynamical axion field. 
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Introduction. — The electrodynamics is a [/(I) gauge the- 
ory usually defined by the Maxwell action with the La- 
grangian density 
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where E and B represent electric and magnetic fields, 
and the permittivity e and permeability fi are both unity 
in vacuum. Gauge invarianee allows an additional term 
in the Lagrangian density 



(2) 



where a is the fine structure constant. Integrating over a 
closed space-time with periodic boundary conditions, we 
obtain the quantization Se = J d'^x Cg = On, where n is 
an integer. Namely, Sg is a topological term. The quan- 
tization also implies that the bulk properties depends on 
9 only modulo 2tt. While Sg generically breaks the parity 
and time-reversal symmetry, both symmetries are intact 
at = and 6 = tt. 

The topological term was originally introduced in high- 
energy physics. In particular, the similar term can be de- 
fined for the quantum chromodynamics, which is a non- 
Abelian gauge theory for the strong interaction. Within 
the Standard Model of Particle Physics, there is a pri- 
ori no reason to set to the time-reversal invariant val- 
ues. If has a generic value, a strong violation of the 
time-reversal or of the CP symmetry should follow, con- 
flicting with current experimental limitations on the CP 
violation. The solution proposed in Ref. [l| to address 
this problem introduces a dynamical pseudo-scalar field, 
called axion field, which couples to E ■ B, so that its 
expectation value absorbs the parameter 6. Let us, for 
the moment, use the same symbol 9 for the axion field. 
It was shown that the axion field would relax into the 
groundstate corresponding to = 0, thereby dynami- 
cally recovering the time-reversal symmetry. This solves 
the "strong CP problem" . At the same time, it leads to 
prediction of a new particle, the axion, corresponding to 
quantum of 0-field. While the axion is a possible compo- 
nent of dark matter, direct detection of the hypothetical 
particle so far remains elusive. 



Condensed matter physics often provides realization of 
intriguing theoretical concepts, which originate in but are 
rather difficult to observe experimentally in high-energy 
physics. The quantum electrodynamics with the topolog- 
ical term ([2]) is such an example. It was pointed out in 
Ref. [H that this system at the nontrivial, time-reversal 
invariant angle = tt is an effective theory for topological 
insulators. In fact, it was pointed out earlier Q that the 
same theory with a smaller value of 9 describes magne- 
toelectric effects in Cr203 [i]. (See also Ref. 0.) 

The topological angle is a static parameter for a topo- 
logical insulator. Nevertheless, in Ref. it was pointed 
out that, when there is an antiferromagnetic order in an 
insulator, the magnetic fiuctuations can couple to elec- 
trons, playing the role of the dynamical axion field. In- 
teresting effects due to the dynamical axion field were 
predicted in the presence of an applied magnetic field. 
Such a system is called "topological magnetic insulator" 
(TMI), and it gives a condensed- matter realization of the 
axion electrodynamics. 

It should be noted, however, that the presence of the 
dynamical axion field does not require the system to be 
a topological insulator; a topologically trivial insulator 
could have the dynamical axion field if there is an appro- 
priate coupling between magnetic fluctuations and elec- 
trons. 

In this paper, we study the instability of the axionic 
electrodynamics in (3-1- 1) dimensions with a massive ax- 
ion fleld, and its possible realization in magnetic systems. 
In particular, we show that the instability leads to a com- 
plete screening of an applied electric field above a critical 
value, proportional to the axion mass. This also leads 
to a spontaneous generation of a magnetic flux from the 
material. We will discuss how this phenomenon can be 
detected experimentally. 

Instability. — We consider the axionic electrodynamics 
defined by the Lagrangian density, in which the electro- 
magnetic field are coupled to the axion field (p, 
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where J , Vi=x.y.x, and m are the stiffness, velocity, and 
mass of the axion The time reversal symmetry is 
broken unless = or tt. For application to magnetic 
materials, the axion velocity Vi can be anisotropic and is 
much smaller than the speed of light in vacuum, which 
is set to unity. 

Suppose we turn on the electric field E in the z- 
direction and consider fluctuations with momentum k in 
the a;-direction. The axion mixes with a photon polarized 
in the y-direction, giving rise to the dispersion, 

UJ^ =i ((C'2 + J.2)fc2 ^ ^2) ^ 



where 



Er: 



m /(27r)3g2j 



a 



(4) 



(5) 



c' = l/y/e/7 is the speed of light in the medium, and 
~ i'i=2- In particular, if i? > i?crit, we find w < for 
the range of momentum. 
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The energy-density H in the material can then be ex- 
pressed in terms of _B as. 
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Minimizing the energy density determines B. 

For example, if we turn off the axion by sending ra^ 
oo, the minimum energy configuration is. 



E = 
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as expected from the Witten effect [7|. We note that, 
although the bulk theory defined by eq. ([3]) has 27r- 
periodicity in 9, the results here are no longer periodic in 



the boundary breaks the periodicity. 
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Namely, the electric field larger than Edit is unstable. 

Consider a fiat plate of the material described by ([3]) 
and sandwich it by (non-topological) insulators with per- 
mittivity eg. Apply a constant external electric field Eq 
perpendicular to the interface between the material and 
the ordinary insulator. The boundary condition at the 
interface gives. 



FIG. 1: The energy density ([T0|) for e = and eo^^o = 
1.3e-Bciit, which is slightly above the critical value eo-Bo = 
eE'crit- The energy 'H(-B) shows a double-well structure and 
the spontaneous breaking of the time reversal symmetry 
B —7- —B,<j> — >■ —(j). The magnetic field corresponding to the 
potential minima is spontaneously generated, resulting in the 
screening of the electric field to iJ = -Bcrit- The magnetic field 
is normalized by the asymptotic value Boo given in eq. (|14|) . 
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where E and B are components of the electric and mag- 
netic fields normal the interface. In addition, the con- 
servation of magnetic fiux enforces that B is continuous 
across the boundary. For now, we assume the Neumann 
boundary condition for the axion field 0, as it enables a 
simple analysis. Later, we will study the case with the 
Dirichlet boundary condition, which is relevant for phys- 
ical realization. 

Assuming homogeneity of the fields, the equation of 
motion sets, 
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The boundary condition ^ then gives, 

_ epEn - a9B/TT 
' ~ l + d^B^/El,,- 
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Let us first analyze the energy density pUj) when 6* = 0. 
When the applied electric field is lower than the critical 
value as in epi^o < e-E-critj the only solution is. 



Eq, B = 0. 



(12) 



However, if the external electric field Eq is raised above 
—Ecrit, we find states with lower energy given by. 



E — i^crit, B — ±-\//i_Ecrit(eo£'o ~ e£^crit)- 



(13) 



Thus, there is a second order phase transition at eoEo = 
eEcrit- Above this value, the electric field inside the ma- 
terial stays constant and the magnetic field is increased 
instead. We show an example of the energy density as a 
function of B above the critical value, in Fig. [1] 

Now, let us turn to the case with 6* 7^ 0. When 
eq-Eo ^ eEciit, there is a unique minimum energy config- 
uration, which reduces to ([TT|) in the limit of Eq 0. For 
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egEo ^ eE'crit, the energy density has two local minima 
in B. Since the time-reversal symmetry is broken explic- 
itly ioT 9 {9 = TT does not preserve the time reversal- 
symmetry in the presence of the boundary) , the two min- 
ima have different energies. In the limit of i?o — >■ oo, the 
more stable of the two behaves as 



E — > i^crit, B — > Boo — \J eoA*^o-Ec: 



(14) 



One can show that the configuration (|TT|) for small 
is smoothly connected to ([Ti|) for large i?o- Namely, the 
phase transition at Eq-Eo = e^-crit is smoothened out when 
6* 7^ 0. However, for realistic values of ^ tt, the smooth- 
ing effect is small and the solution is similar to the one at 
= with a second-order phase transition, except that 
the groundstate is chosen uniquely for any £^o- 
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FIG. 2: The electric field E (red) inside the material and 
the induced magnetic field B (blue) as functions of the ex- 
ternal field i?o, evaluated at S = vr. B and E are deter- 
mined by the global minimum of the energy density (|10|l . 
The fields i5o , E and B are normalized by ^ -Bcrit , Scrit , and 
Boo(Eq = 2^i?crit), respectively. The behavior appears to 
indicate a second-order phase transition, although the transi- 
tion is smoothed out by the effect of non-zero Q. 



We note that, in realistic systems, higher order terms 
of are expected in the effective theory, since the range 
of ordered magnetic moment, which corresponds to shift 
of (/), is bounded. However, although it will be impor- 
tant for E ^ -Ecrit, this does not affect the existence of 
the transition, because (/), i? ^ in the vicinity of the 
transition. 

For experimental realization of the screening effect dis- 
cussed below, the Dirichlet boundary condition <^ ~ 
will also turn out to be relevant. In this case, the solu- 
tion is no longer uniform in the z-direction. Nevertheless, 
away from the boundary, the solution asymptotically ap- 
proaches to the stationary configuration, which is given 
by the solution obtained in the above for the Neumann 
boundary condition, with the replacement 9 — > — 0o- The 
screening of the electric field occurs in the transient re- 
gion with a finite length. 

If we ignore the coupling to the electromagnetic field, 
the axion field approaches to the stationary solution 



exponentially, with the decay length of vjm. However, 
we should take into account mixing of the axion and the 
photon via the coupling (pEB. This gives the axion field 
the effective mass moff, which is given by 
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When <j>o = —9 = and the applied electric field is above 
the critical value eg -Bo > eBcrit, the magnetic field is 
given by ([T5|) . When 0o ^ 0, this is the asymptotic value 
for Eq — oo. Substituting this value of B in the above, 
we find that the screening occurs within the lengthscale 
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Physically, the screening of the electric field occurs be- 
cause of the induced charge density oc V0- B in the axion 
electrodynamics Namely, the axion field </> is shifted 
inside the material creating the gradient V0 near the 
boundary. By generating the magnetic field B, this in- 
duces a charge density at the boundary, screening the 
electric field inside the material. 

From this physical picture based on screening, we can 
also understand why the second-order transition occurs 
when (f>o = 0. Under a given applied electric field, the 
sign of the charge needed to screen the electric field is 
uniquely determined. However, induced charge is pro- 
portional to dz<f> ■ B, which is invariant under the time 
reversal </)—)■ —(t>,B — > —B. The symmetry is preserved 
for the boundary condition 0o = 0, and there is no pre- 
ferred sign of B. The system breaks the symmetry spon- 
taneously to produce screening charge, for E > Bcrit- 

On the other hand, (j)Q ^ introduces a gradient of cf) 
near the boundary, choosing the preferred sign of B to 
produce screening charge. Thus the symmetry is broken 
explicitly and the phase transition is smeared. 

Possible Experimental Realization in Magnetic 
Materials. — Let us discuss possible realization of 
the axionic instability in condensed matter systems. 
In Ref. Bi2Se3 doped with 3d transition metal 
elements such as Fe (Bi2Se3-Fe hereafter) is discussed as 
a candidate for TMI, which is described by the axionic 
electrodynamics. Although the mechanism proposed in 
this paper is not restricted to any particular system, 
we shall examine possible physical realization using 
Bi2Se3-Fe as a reference. Because of the magnetic dop- 
ing, a magnetic order M~ , which is ferromagnetic in the 
xy-plane and antiferromagnetic along the z-direction, 
may appear Q. In the following, we assume that the 
electric field is applied along the z-axis. 

The (relative) permittivity, axion mass, and axion cou- 
pling in Bi2Se3-Fe are estimated 0] to be e ~ 100, 
?7i - 2 meV, and = a(0.4 T/meV)^. This 

yields a rather high value of Bcrit = 2.4 x 10^ V/m, which 
is above the breakdown field of typical semiconductors. 
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The critical field i^crit is reduced for smaller axion mass 
m. In the ordered phase, the axion mass is proportional 
to the spontaneous magnetic order M~ 0. Thus, the ax- 
ion mass may be reduced by tuning the system close to 
the critical point, so that the antiferromagnetic order be- 
comes small. If, for example, i^crit ~ 10^ V/m and a TMI 
film of thickness ^ 10 ~® m is used, the voltage difference 
across the system is of the order of 0.1 V. The corre- 
sponding energy is well below the band-gap ^ 0.3 eV of 
Bi2Se3, justifying the low-energy description. 

However, TMIs such as Bi2Se3-Fe arc expected to have 
surface electronic states. For the undopcd topological in- 
sulator (such as Bi2Se3), the surface states are described 
by massless Dirac fermions. In the doped TMI, the sur- 
face state has a gap due to the time-reversal sym- 
metry breaking. In Bi2Se3-Fe, it was estimated 0] that 
ms = 1 meV and m = 2 meV; these are of the same or- 
der of the magnitude. However, in order to suppress the 
screening effect by the surface states and to enhance the 
effect by the dynamical axion in the bulk, it is desirable 
to keep the surface Dirac mass ms large while tuning the 
system close to the criticality to reduce the axion mass 
m. This could be achieved by sandwiching the TMI by 
ferromagnets, which enforces the magnetic order near the 
surfaces and determines the boundary value 4n) . 

The axion mass m and the spin-wave velocity v are 
expected to be of the order of U and Je^-a, respectively, 
where Jex is the effective exchange interaction and a is 
the lattice constant. Thus, the screening length is 
expected to be of the order of the lattice constant when 
U is sufficiently large. This implies that the present effect 
could be observed in thin film samples. 

Experimentally, the instability discussed in this paper 
is observed as an increase of capacitance above the criti- 
cal electric field. Another experimentally observable con- 
sequence is the generation of magnetic field as in eq. (fT4| . 

In the above scenario, the necessity to gap out the sur- 
face Dirac mode by sandwiching the TMI by ferromag- 
nets may pose an additional complication. If we begin 
with a topologically trivial insulator, there is no surface 
Dirac mode. In fact, for the instability discussed here, it 
is not essential that the system is based on a nontrivial 
topological insulator. Even in a topologically trivial in- 
sulator, if there is an appropriate coupling with the mag- 
netic order and the electrons, the magnetic fluctuations 
can play the role of the axion field, following the argu- 
ment parallel to that of Ref. @. The difference is that, 
starting from a trivial insulator we expect 0^0, instead 
of 6* ~ TT in a TMI. However, as long as the dynamical 
axion <j) is present in the low-energy effective theory ([3|. 
the mechanism proposed in this paper should work. The 
class of magnetic materials with a dynamical axion field 
may be called axionic insulators. 

Cr203 is a topologically trivial insulator with the band 
gap of 3.4 eV, and a magnetic long-range order. As 
we have mentioned, it exhibits magneto-electric effects. 



which can be described by the axion electrodynamics. It 
is also suggested that the antiferromagnetic fiuctuations 
in Cr203 play the role of the dynamical axion field 0]. 
The present mechanism could in principle be realized in 
such a system. The large band gap itself is advantageous 
for observation of the effect, as the material can with- 
stand stronger electric field. However, the large band gap 
also implies a large coupling g, which is approximately 
proportional to the band gap. This results in a larger 
critical field because of eq. ((S]), cancelling the advantage. 

The axion mass m is proportional to the on-site 
Coulomb repulsion U and the magnetic order The 
Coulomb repulsion U in Cr203 is estimated to be about 
5 eV [l3|; thus the axion mass m is expected to be of 
the same order, resulting in a too large value of i?crit- In 
order to observe the instability, we need to find a system 
with smaller axion mass. Again, this could be achieved 
by tuning the system near the magnetic criticality, for 
example by mixing with nonmagnetic ions. 

We would like to end with a remark on earlier papers, 
where similar instabilities were discussed. It was pointed 
out in Ref. 11 , [l^l that the Maxwell theory with the 
Chern-Simons term in (4+1) space-time dimensions is 
unstable in the presence of a constant electric field. Sub- 
sequently, a similar instability was also found in its di- 
mensional reduction to (3 -|- 1) dimensions: the massless 
axion electrodynamics [l3,[l3|- (See also the earlier work 
[ist.) In this paper, we extended these results for the 
(3 -f- 1) dimensional theory with a massive axion field, 
clarified the eventual fate of the unstable gauge theory 
in the Minkowski space, and proposed its possible exper- 
imental realization. Although these earlier papers were 
in the context of the AdS/CFT correspondence, nothing 
in this paper assumes the AdS/CFT correspondence or 
depends on results obtained by it. 
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